Some integrals related to the Fermi function 



o 

(N 



in 

(N 



43 

^ : 

-i— > 
o3 



> 

O 

(N 

rn 
O 
m 



G. Bevilacqua 

Department of Information Engineering and Mathematical Sciences, 
University of Siena, Via Roma 56, 1-53100 Siena, 



Some elaborations regarding the Hilbert and Fourier transforms of Fermi function are presented. 
The main result shows that the Hilbert transform of the difference of two Fermi functions has an 
analytical expression in terms of the ^ (digamma) function, while its Fourier transform is expressed 
by mean of elementary functions. Moreover an integral involving the product of the difference 
of two Fermi functions with its Hilbert transform is evaluated analytically. These findings are of 
fundamental importance in discussing the transport properties of electronic systems. 



I. INTRODUCTION 

In many problem in electronic transport calculation 
the Hilbert transform (HT) of the difference of two Fermi 
function (FF) is required^. Usually this task is accom- 
plished in a numerical way using the Matsubara expan- 
sion of FF or other approximate expansions with faster 
convergence behaviour— ~—. While these methods as well 
as the Sommerfeld expansion are of valuable help in nu- 
merical work, the importance of exact results can be 
hardy overlooked. 

In this paper I will show that the Hilbert and Fourier 
transform of the difference of two FF can be expressed in 
a closed analytical form by means of the digamma func- 
tion. Moreover an universal integral needed in the theory 
of electronic transport in presence of phonon scattering 1 
is exactly evaluated. 

The paper is organized as follow: in Scction|TT]a quickly 
description of the contour integration details used in the 
work is given; in Section [TTT1 the main results are derived, 
while in Section IIVI an universal integral needed in the 
discussion of the asymmetric contribution to the current 
flowing in a molecular wire is worked out. Section [V] 
contains the conclusions. 



II. HILBERT TRANSFORM 



Given a function g(x) its HT is defined as 



H(g)(y) = g H (y) 



-P 



+ OC 



g(g) 
y-x 



dx 



(1) 



where P is the usual Cauchy principal value. 

One powerful technique often used in the evaluation 
of such kind of integrals extends the function involved in 
the complex plane and makes use the residua theorem. 
The typical situation is shown in Fig. [1] 

By the residua theorem the integral over the closed 
circuit, done in counter-clockwise way is proportional to 
the sum of the internal residua, while the integral over 
the real axis interval [A, B] U [C, D] in the limit of R — > 
oo, and, separately r — > gives the HT desired. More 




Figure 1: Circuit in the complex plane used to evaluate 
Hilbert transform. The x symbols denote the poles of the 
integrand located at x n points. 



precisely 



L 



full circuit 



2ni \_. Res(x„) 



r J[A,B]U[C,D] Jj 



(2) 



If the integral over the T half-circumference vanish in 
the limit R — > oo one has 



2-Ki ^2 Res(x„) = J 



[-oo,B]U[C,+oo] J-y 



(3) 



while in the limit r — >■ the first integral in the right- 
hand-side of ([31) becomes essentially the principal value 
and the latter can be related to the residuum in y, i.e. 



-Tii Res(y), where the minus sign is due to the 



clockwise orientation of the path. 

In summary one arrives at the following result 



-P 

^ ./-oo y~ x 



Ax = 2 i Res(x„) + i Res(y) (4) 



where x n are the poles of g(x) in the 3(x) > half-plane. 
In case g{x) has no poles on the real axis the preceding 
formula simplifies as 

1 P f +OC l&L dx = 2iY, Res(x„) - i g{y). (5) 



7T 



y-x 



2 



III. HILBERT AND FOURIER TRANSFORM 
OF THE DIFFERENCE OF TWO FERMI 
FUNCTIONS 

Let apply the above machinery to the difference of FF 
1 1 



fi(x)-f 2 (x) (6) 



where the chemical potentials /ii and inverse tempera- 
tures Pi can be thought as parameters. The poles struc- 
ture of g is composed of two infinite series i) x n — 
Hi + (n + 1/2) 2m/ fix, n = 0, ±1, ±2, . . . with residua 
-l//3i; ii) x m = M2 + (m+l/2) 2m/p 2l m = 0, ±1, ±2, . . . 
with residua 1 //?2- Only those with n,m > are re- 
quired. Collecting all the contributions one obtains 



fJH 



(y) = +2i J2 

oo 

+ 2* £ 



[X X n ) \X X n ) 



y-x 



y-x 



m=0 

i[fi(y)-f2(y)}. 



y-x 



y-x 



(7) 



Noticing that only the first and the fourth term in the 
square brackets in the right-hand side are different from 
zero and equal respectively to (—l/Pi)/(y — x n ) and 
(— l/p 2 )/(y — %m), after some straightforward manipu- 
lations one obtains 



9h 



-I/ft 



y-Mi-(n + l/2)27ri/j8i 

-I//82 

y - - (n + l/2)2m/fo. 

i[h{y)-h(y)]. 



(8) 



Introducing the quantities lOj = 1/2 + (y — fj,j)/2n, 
the sum can be recast in a form that shows the relation 
with the Euler digamma function 'J 6, ' 



E 

71 = 



1 



1 



n + W\ n + w 2 



= *(io 2 )-*(wi) (9) 
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Figure 2: (Top) Difference of two Fermi functions. (Bottom) 
Hilbert transform (solid line) and first derivative (dashed line) 
of the top panel. The parameters are /11 — 0, fi2 — —15, 
01 = 6 and (3 2 = 0.4. 



Various properties of gn can be inferred by inspection. 
For instance 

gH (y) = l\n^ + ^—^ + 0(l/y 2 ) M^oo(12) 
7T Pi y 



Moreover if Pi = p 2 , g(x) is even with respect to the 
point —(fJ.i + /i 2 )/2 while gn{y) is odd. Finally note that 
in the low temperature limit p — > 00 from the asymptotic 
expression of the digamma function 6 



The term fi (y) — f 2 (y) can be further elaborated exploit- 
ing the relation between the ^S> and the FF 

m = l + -S[*(l/2-iz/2n)} (10) 

so after some manipulations and remembering that 
*(z*) = iff(z)* one finds the analytical result 

g H {y, mi,m 2 ,/?i, AO = -K [* (w 2 ) - * (wi)} . (11) 



In Fig. [3 a particular case is shown. 



tf(z) = ln(z)-_ + 0(z- 2 ) 



the elementary result 



9h(v) = - In 

7T 



2/ - M2 



y - Mi 



+ o(/r 2 ) 



(13) 



(14) 



is obtained. 8 

The Fourier transform of g(x) can be easily computed 
with a slight modification of the above procedure. In fact 
choosing A > and using the Jordan lemma closing the 



3 



contour in Ss(x) > one obtains 



d) x n = V + (n + l/2)2m/0, n = 0, 1,2, . . . with 
residua 



9f(X) 



dxe tXx g{x) 



1 p i A u 2 



1 p i A /ii 



sinh(7rA//32) sinh(7rA//3i 



(15) 



r v (n) 



_ z n+l/2 x 



The same expression is obtained choosing A < and 
closing the contour in < 0. Notice that from this 

result one can easily evaluate the Fourier transform of 
gu by mean of the general relation 8 



27T/3 

*(« i>_ - n) + v_ + n + 1) 



t> + - n) - i w+ + n + 1) 



(22) 



+ 00 



The integral is thus expressed as 

+ OC 



dxe lXx g H (x) = -i sign(A) g F (A) . (16) #0) = 2™ ^ [r +w (n) + r_ w (n) + r (n) + r y (n)] 



IV. EVALUATION OF AN INTEGRAL 

Let consider the special case gn(x;u!, —u>, j3, 0) and 
g(x; 0, —V, /3, /3), the integral 



I = 



+ 00 



g(x; 0, -V, 0, p) g H (x; u, -u, /3, 0)dx, (17) 



describes the asymmetrical contribution to the current 
flowing in a molecular wireA. To avoid convergence prob- 
lems it is better to evaluate the Fourier integral 



K(X) 



+ OC 



e lXx g(x; 0, -V, f3, 0) g H {x;oj, -u, (3, 0)dx, 

(18) 

from which one obtains I = K(X — > 0). In fact again 
using the Jordan lemma with A > one needs the poles 
structure of the integrand in > 0. Let us define 

some intermediate quantities to simplify the notation 
c = coth(7ru ), c± = coth(7ri>±), ri± LU = -(±c + c T )/2 
where vq — /3u>/2tt and v± = f3(V ±w) /2tt. Moreover let 
z = cxp(— 2irX/f3) < 1. The poles are found straightfor- 
wardly: 

a) x n = uj+(n+l/2)2iri/ f3, n = 0, 1,2, . . . with residua 

r +u (n) = ^e iXu z n+1 / 2 r, +u (19) 

b) x n = -oj + (n + l/2)2m/P, n = 0, 1, 2, . . . with 
residua 



(20) 



c) x n = (n + l/2)2m/0, n = 0, 1, 2, . . ., with residua 



r (n) = — 2 n+1/2 3 [*(1 + n + i u ) + « «o - n)l 
7rp 



(21) 



(23) 

The first two terms are trivial geometric series whose 
result is 



1 



n=0 

i 



(24) 



2/3 sinh(7rA//3) 



fa 



For the third term using the results of Appendix [A] one 
gets 



+ 00 



1-2 



Ro = ^2r (n)= lZ i$S[S+{ivQ,z)] 

n=0 ^ 

similarly for the last term one finds 



(25) 



+ 00 



e av z i/2 
R v = }^ rv(n) = — x 

S+(-iv-,z) - S + (-iv +1 z) + i^- — — 

2 1 — z 

(26) 

The integral is now completely solved as 

A' (A) = 2vr i (R u + R + R v ) . (27) 

To obtain the limiting value A^(0) one needs to care- 
fully extract the divergent terms from R u , Rq and Ry 



Ru = 



2tt 



A 



Rq + R\ 



2ir 

i (r] +UJ - ?7_^) 1 
2tt A 



2h(v ) + h(v-)-h(v+) 



— (c_+ c+ ) + 0(X), 



(28) 



(29) 



4 



where the property 6 - ^(^(ix)) = 1/(2 x) + 
(tt/2) coth(7ra:) is used and the function h(x) is 

h{x) =xU[^{ix)]. (30) 

As expected the divergent terms cancel each other and 
one gets the hnal result 



I 



h(v + ) - h(v-) - 2h(v ) 



(31) 



Using the asymptotic expansion 6 - h(x) = x hn(x) + 
l/(12x) + l/(120a; 3 ) + 0(x~ 5 ) the limit of low tempera- 
ture is obtained 



7=1 

7T 



(V + u) ln(V +u)-(V-w) ln(V -u))-2u> In lj 
2ttV 2 



u (V 2 - u 2 ) 



p-' + o(p-*) 



(32) 



where the first term is related to the elementary case 
(j!4|). On the other hand for small x one has^ h(x) — 
*(1) + C(3) + 0{x 5 ) (C(3) w 1.202056) and the limit 
of high temperature follows 



12 C (3)a,y3 

(2tt) 3 P lP J 



(33) 



The double sum can be decoupled using the rule 



-f-oo n 



-\-oo + 00 



(A4) 



thus obtaining 



S+Kz) = — [*H (A5) 

where the last term is a particular case of the Lerch 
transcendent! cf)(z,w) — $(z, 1, w) 



$(z,s,w) = E 



— ' (n + w) s 



(A6) 



Notice that <j> can be also expressed in term of Gauss 
hypergeometric function^ 



(z,w) = — 2 F 1 (l,w:l + w;z), 



(A7) 



and has a logarithmic singularity in z = 1 as can checked 
from the expansion 6 



V. CONCLUSIONS 

The paper derives analytical results for some inte- 
grals related to the difference of two Fermi functions: its 
Hilbert transform is evaluated in term of the digamma 
function while the Fourier transform involves only ele- 
mentary functions. Additionally an exact result useful in 
the discussion of the asymmetric part of the electronic 
current in presence of phonon scattering is reported. 



Appendix A: A summation formula 

The series 



6(z, w) = V (1-2)" [*(n + 1) - 9(w + n) - ln(l - z) 

n=0 

which is valid if |1 — z\ < 1 and | arg(l — z)\ < it. 

From the reflection property of the digamma \I/(1— z) = 
*(z) + 7rcot(7rz) follows that *(w-n) = *(w + n+l)- 
7rcot(7rw) and thus 



*(w - n) z n = S + (-w,z) - 



7T COt(7T Vj) 



(A9) 



n=0 



Moreover 



S+(w,z) = ^*(w; + n + l)z" |z|<l, (Al) J2W(-w-n) + V(w + n+l)]z n = 2S+(w,z)+tt^ 



C0t(7T w) 



n=0 



can be expressed in a closed fornA In fact exploiting the 
telescoping properties of Of 

" 1 

*(w + n + l)=V- + m(w) (A2) 



k + W 

k=0 

and multiplying by z 11 and summing one gets 

+oo n n 

S + (w,z) = - *(«,) + EEirh- ( A3 ) 

1 — z z — ' k + w 

n=0 fe=0 



— z 
(A10) 

When z — > 1~ along the real axis, the behavior of 
S+(w, z) is easily extracted from (|A8[) 



\I>(1) - ln(l - z) r , , 
S+(w, z) = K z '- + w [#(2) - $(w + 1)] 



1 - z 

wln(l-z)+0(l-z). 



(All) 
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